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ABSTRACT 

We presenl a numerically stable and computational simple 
method for fast and robust measurement of optical channel 
gain. By transmitting adaptivcly designed signals through the 
channel, good accuracy is possible even in severe noise condi- 
tions. 



1. INTRODUCTION 

The measuring of opUcal channel gains Is a key element in 
many applications. Measuring channel gains means determin- 
ing the change in intensity when a signal is transmitted from 
an emitter to a feceiver. A well-known and simple applica- 
tion is an automatic door, which responds whenever a person 
Is reflecting the emined signal, and thereby significanily in- 
creasing the channel gain. Another example is measuring the 
thickness of paper. A more subtle example is determination 
of spatial position by comparing the intensities of a multitude 
of reflections from a single object. A typical way of making 
this type of measurements is emitting a simple signal, such as 
a harmonic or square wave signal^ since they are both easily 
constrvcied and measured with analog electmnics. Such solu- 
tions have two major disadvantages: The signals are sensitive 
to frequency located disturbances, and it is difficult to detect 
and avoid/neutralize such disturbances; ^ 

We propose a measuring method which is highly accurate 
in moderate noise condiUons. and less accurate, but very ro- 
bust, in severe noise conditions. This is achieved by using 
two closely related digital signal design algorithms; a 'Tjcst 
case" and a "worst case" algorithm. The former is based on 
the wavelet transform (WT). while the latter is based on the 
Rudin-Shapiro transform (RST). They arc both simple, numer- 
ically stable, and post-processing friendly making them ideal 
for implementation e.g. in a fixed point DSP. By intixxiucing a 
signal processor it becomes possible to continuously redesign 
the signals for improved SNR, and thereby maintaining the ac- 
curacy in changing and/or severe noise conditions. 

2. DESIGNING THE DIGITAL SIGNALS 

The two design algorithms are based on the wavelet packet 
transform scheme; it is fast, numerically stable, works well 
in fixed point arithmetics, and has low program complexity. 

"niis woric is in part iupportcd by the Danish Technical Science 
Foundation Grant no. 9701481 . 
Patent pending. 



The best case algorithm uses the classical WT to create sig- 
nals which are near-orthogonal to expected noise occurrences, 
while the other algorithm uses the RST to create an all-spectrum 
signal, which by nature has low sensitivity with respect to time 
and frequency located noise occurrences. The diflTercnce is 
essenUally that the WT algorithm "searches for holes" in the 
current noise, while the RST algorithm spreads information in 
time and frequency to reduce the impact of localized distur- 
bances. The prefered method depends on the noise conditions. 
If there are casy-to-find holes in the noise, the former can pro- 
vide very accurate measurements. If, however, the noise is dif- 
ficult to define or is changing rapidly, the latter method pro- 
vides lesser accurate, but more robust measurements. 

A good introduction to the wavelet theory is Wickerhauser 
[6]. A mathematically rigorous treatment of the subject is given 
In Daubcchies [3]. For more material on Rudin-Shapiro poly- 
nomials see Brillhan (1 ] 

2.1. The wavelet transform 

The WT based algorithm takes a simple, time localized sig- 
nal (see figure I for an example), and inversely wavelet packet 
transform it, which results in a frequency localized signal. A 
typical maximum spread is also shown in figure 1 . After trans- 
mission the signal is forwardly transformed to reproduce the 
original, now noisy simple signal, and by inner prx>ducl with 
the "clean" original signal the u^nsmission intensity (the chan- 
nel gain) is determined. Since the original signal is completely 
known, it is also possible to obtain an estimated accuracy of 
the channel gain measurement. This is accomplished by taking 
inner product between the transmitted, transformed signal and 
a number of signals orthogonal to the original signal. If these 
quantities are small the transmission was most likely subject to 
only mild noise. This trick provides an easily calculated guide- 
line to how much one can trust the current measurement. If 
each measurement is vital a number of signal restoration pro- 
cedures (not further described here) can be applied. These also 
benefit from the complete knowledge of the original signal. 
Note that the small number of non- vanishing coefficients of 
the original signal in all cases significantly reduce the amount 
of calculations. 

Because the transform is linear and has perfect reconstruc- 
tion, it is also easy to make a good sample by sample estimation 
of the noise, which provides valuable information on any dis- 
turbances. This makes it possible to do real time adaptation 
of the signal. The mchlods combines the ability of the WT 
to produce predefined trade-offs between time and frequency 
information with the freedom in design of the original signal. 
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Figure I: Uppermost a simple signal (here a sampled chirp) 
Below (in solid) the absolute value of the Fourier iransform 
of the 3-scaIe inverse wavelet packet transform with symlets 
6 [4], The dashed curve shows the maxima! frequency spread- 
ing for any (suitably normalized) signal with coefficients van- 
ishing outside [128; 159], 



Thereby it is possible to adapt the method to virtually any type 
of noise, in particular disturbances with large temporal extent 

The Rudin-Shapiro transform 

The RST is defined through a slighUy extended version of the 
remaricable Rudin-Shapiro polynomials, introduced in 1951 by 
H. S. Shapiro in his master's thesis, and published in 1959 by 
Rudin [5]. Define the polynimals 



(I) 



with Sm € {0, 1}. It immediate follows that for all \z\ ~ 1 

Consequently, 

max |/>m(e'<)| < v/2||P^(e'<)||2, (2) 

guaranteeing a certain flatness of the polynomials. A construc- 
uon similar to (1) is found in Byrnes [2]. The coefficients 
of the polynomials can also be constructed with the Rudin- 
Shapiro transform, which is really a modified wavelet packet 
Haar transform. Define the unitary transform H„ : R^* 
R . n > 1, as 

f f* 1 = -L (-1)* 1 r 1 

bk+J-'J v/2 ll -(-1)*J 

for * = 0, . . . , 2"-' - 1 when mapping x to y. Then, wiOi 
lin being the inverse of H„, 



H„ 
0 



0 



H„ 



transforms the canonical basis for into the coefficients of 
the 2 possible Pa'(*)'s. Hence these coefficients consUiuies 
an onhonormal basis, which not only consists of only ±t, but 
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Figure 2: Uppermost is the result of applying the inverse 
Rudin-Shapiro transform to the 29th canonical basis vector of 
length 64. Below the frequency response (the absolute value of 
the cofTCsponding P^{e*^) polynomial), 

also, due to (2), has a remarkable frequency response. In fig- 
ure 2 is an example of such a basis clement and its frequency 
response. The RST based algorithm is applied in much the 
same way as the WT method. A simple signal is inversely 
transformed prior to transmission. This will produce an all- 
spectrum signal. Upon transmission the signal is forwardly 
transformed yielding the original, simple signal with noise. 
The post-processing is equivalent to that of the WT method. 

3. CONCLUSION 

TWO computationally simple and numerically robust algorithms 
for measuring optical channel gains were presented. One of 
the algorithms provides excellent accuracy in moderate noise 
conditions, while the other has reduced, but very robust, ac- 
curacy even in severe noise condiUons. Combining the two 
algorithms, cither by applying the most suitable one. or jointly 
in two parallel systems, is easy due to their similar program 
and computational structure, and the result is a versatile op- 
tical channel gain method. The low complexity and numer- 
ical stability of the wavelet packet transform scheme and of 
the post-processing (mainly inner products) also makes this ap- 
proach fast and suitable for low cost hardware implemcntaUon. 
For the commercial aspects of these methods, please refer to 
www . beaxncon t ro 1 . com. 
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